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ABSTRACT 

On December 27, 2004, a giant flare from SGR 1806-20 was detected on earth. Its thermal spectrum and 
temperature suggest that the flare resulted from an energy release of about 10 47 erg/sec close to the surface of 
a neutron star in the form of radiation and/or pairs. This plasma expanded under its own pressure producing a 
fireball and the observed gamma-rays escaped once the fireball became optically thin. The giant flare was followed 
by a bright radio afterglow, with an observable extended size, implying an energetic relativistic outflow. We revisit 
here the evolution of relativistic fireballs and we calculate the Lorentz factor and energy remaining in relativistic 
outflow once the radiation escapes. We show that pairs that arise naturally in a pure pairs-radiation fireball do not 
carry enough energy to account for the observed afterglow. We consider various alternatives and we show that if 
the relativistic outflow that causes the afterglow is related directly to the prompt flare, then the initial fireball must 
be loaded by baryons or Poynting flux. While we focus on parameters applicable to the giant flare and the radio 
afterglow of SGR 1806-20 the calculations presented here might be also applicable to GRBs. 
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via annihilation of the magnetic field of a highly magnetized 
neutron star. This annihilation deposits energy at the form of 
photons and pairs near the surface of the neutron star. The pair- 
radiation plasma evolves as an accelerating fireball (Goodman 
1986; Paczynski 1986; Shemi & Piran 1990; Piran et al. 1993; 
Meszaros et al. 1993; Grimsrud & Wasserman 1998) resulting 
in a thermal radiation burst carrying the bulk of the initial en- 
ergy with roughly the original temperature and a fraction of the 
energy in the form of relativistic pairs. The thermal spectrum 
of the flare and its temperature support this picture. 

Here we compare the energy required to produce the radio 
afterglow (calculated in §2) with a new simple calculation (§3) 
of the energy of the pairs outflow. As the available pairs energy 
is short by at least two orders of magnitude we consider (§4) 
baryonic or electromagnetically loaded fireballs, again provid- 
ing new simple estimates for these cases. We compare our 
calculations of the fireball evolution to previous works in §5 
and find that our simple estimates correct previous works. Fi- 
nally, we discuss the implications of our results to the giant flare 
of Dec 27 th flare from SGR 1806-20 and of August 27"' from 
SGR 1900+14 in §6. 

2. A LOWER LIMIT ON THE RADIO AFTERGLOW'S ENERGY 

The energy emitted in the radio during the second week af- 
ter the burst is about 10 38 erg (Cameron et al. 2005; Gaensler 
et al. 2005). This is clearly an absolute lower limit to the en- 
ergy of the relativistic ejecta. A stronger limit can be obtained 
if we assume synchrotron emission as suggested by the op- 
tically thin spectrum, F v oc v~° 1 and the observed linear po- 
larization (Gaensler et al. 2005). We employ here the famil- 
iar equipartition method (Pacholczyk 1970; Scott & Readhead 
1977). To emphasize the robustness of this method, which 
depends only on the assumption of synchrotron radiation, we 
sketch the derivation here. We characterize the emitting region 
by the number of its pairs N, the magnetic field, B, and the typi- 
cal thermal Lorentz factor of the pairs, j e . We ignore the mildly 

1 Through out the paper we consider a distance of 15kpc even though SGR 1806-20 might be closer by a factor of <2 (Figer et al. (2004); Nakar et al. (2005); 
Cameron et al. (2005) see however McClure-Griffiths & Gaensler (2005)). We express the dependence of our calculations on the flare energy as implied from this 
distance 

2 A radio afterglow was detected also after the August giant flare from SGR 1900+14 (Frail et al. 1999) 
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1. INTRODUCTION 

The giant flare from SGR 1806-20 was the most powerful 
flare of gamma-rays ever measured on earth. It lasted about 
0.2sec. Its fluence of w 1 erg cm" 2 (Hurley et al. 2005; Palmer 
et al. 2005), correspond to energy of 3 x 10 46 erg released at a 
distance of 1 . Our conclusions do not depend on the exact dis- 
tance. 15kpc (Corbel & Eikenberry 2004). This energy exceeds 
the energies of the giant flares from SGR 0526-66 (the famous 
March 5"' 1979 event) and from SGR 1900+14 (August 27"' 
1998) by a factor of a hundred (Mazets et al. 1999). The spec- 
trum of the flare is consistent with that of a cooling blackbody 
spectrum with an average temperature of 175 ±25keV (Hur- 
ley et al. 2005). The reflection from the moon, as detected by 
Helicon-Corona-F (Golenetskii et al. 2004) with a fluence of 
7.5 x 10~ 7 erg cm" 2 , provides an alternative mean of estimating 
the fluence at the energy range 25-400keV. The albedo of the 
moon in this energy range is around 0.25 (Nakar et al. 2005), 
resulting in an isotropic energy release of about 10 46 erg. This 
value is consistent with the fluence and the spectrum measured 
by RHESSI, since for a 170 keV black body, only a quarter of 
the energy is radiated in the 25-400 keV range. 

Like the two other giant flares this flare was followed by 
a pulsed softer X-ray emission that lasted more than 380sec 
(Mazets et al. 2004). Radio afterglow was detected from VLA 
observations (Cameron et al. 2005; Gaensler et al. 2005) 2 . Af- 
ter one week the radio source was extended with a size of 
9 = 40-80masec corresponding to a radius of 0.6-0.9 x 10 16 cm 
and an average velocity of 0.3c-0.5c (at d = 15kpc). Therefore 
a significant amount of energy was emitted in the form of a rel- 
ativistic ejecta around the same time that the 7-rays flare was 
emitted. 

According to the standard magnetar model (Duncan & 
Thompson 1992; Paczynski 1992; Katz 1994; Thompson & 
Duncan 1995, 1996; Woods et al. 2001; Eichler 2002; see also 
Katz 1982 for related early ideas) , the giant flare is produced 
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relativistic motion. The most conservative assumption that all 
the electrons emit at a synchrotron frequency of v R = 8.5GHz 
requires that the pairs have a Lorentz factor of: 

While to obtain the observed flux the number of pairs must sat- 
isfy: 

a T m e c z B 

The sum of the energies of the pairs and of the magnetic field, 
E = R 3 B 2 1 6 + Nm e c 2 j e , is minimized once the former is 3/4 of 
the latter: 



by the condition that their annihilation rate n±a(f3±)f3±c equals 
the local expansion rate c/(R/j) w c/Rq: 



N± = 



4nR ctT^ 

<7 T Tl 



<2Jxl<t*E%%?tf, (6) 
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We assumed here, conservatively, that all the emission is radi- 
ated at vr. As the observed afterglow shows a spectral power- 
law over more than an order of magnitude in frequency our es- 
timate for E min should be larger by a factor of 5 (Pacholczyk 
1970; Cameron et al. 2005; Gaensler et al. 2005) bringing it to 
- 3 x 10 43 erg. 

If the radio afterglow of this SGR arises, like in GRBs or 
SNRs, from shocks going into the surrounding medium, farther 
constraint exists. The observed shock size determines its ve- 
locity and therefore the typical electron Lorentz factor. The ex- 
ternal density and the Lorentz factor dictate the magnetic field. 
These additional constraints will force j e and N to deviate from 
equations 1 & 2 and therefore result in a higher energy estimate. 

3. PAIRS-RADIATION FIREBALL 

Consider an energy E = 1 46 erg that is deposited in the vicin- p . 

ity of a neutron star (R « 10 6 cm) as photons with a typical en- 
ergy e 1 k> 500 keV (with a black body temperature of 170 keV, . . . . 3 
more than half of the energy is in photons with e > e 7 ). If and their klnetlc ener gy ls " : 



where n± is the pairs density and T±, the temperature in the fire- 
ball rest frame at R±, is ~ 18keV (Goodman 1986; Paczynski 
1986; Shemi & Piran 1990). We have used here the fact that at 
low energy the cross section for annihilation is a((3±) ~<jt / '/3± 
with (3±c = v± the thermal velocity of the pairs in the local 
frame. The photons decouple from the pairs around the same 
time that the pairs freeze out. While at this stage the thermal 
velocity of the pairs v± is much smaller than c, the cross sec- 
tion for annihilation is larger than the cross section for pho- 
ton scattering by the same factor. The escaping photons, which 
carry most of the energy, have a quasi-thermal spectrum (Good- 
man 1986; Grimsrud & Wasserman 1998) with T b s = TT± = To, 
roughly the initial temperature. 

Even after the pairs stop annihilating and the photons de- 
couple from the pair the pairs do not decouple from the pho- 
tons. The huge photon flux continues to accelerate the pairs 
(Meszaros et al. 1993; Grimsrud & Wasserman 1998). The ac- 
celeration continues as long as the force that the photon field 
applies on an electron (positron) is sufficient to accelerate the 
electron so that it remains at the same Lorentz factor as the bulk 
of the photon field. The condition for effective acceleration is 
that during the time that the radius doubles the work done by the 
photon field on an electron Eaj /AirRct^ 2 is larger than 7m c c 2 , 
the energy that the electron needs to gain during this period in 
order to keep up with the accelerating flow. This implies that 
the pairs accelerate until they reach a bulk Lorentz factor of : 
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the energy is released instantaneously then the duration of the 
observed emission would be Rq/c w 0.1ms/?o, 6 - To be com- 
patible with the observed duration of the flare, the source must 
have been active for a time comparable to the observed dura- 
tion, t w 0.1 sec. The optical depth for pair production would 
be: 



E± =N±m e c z T ± « 1.4 x 10 41 erg E 46 R Q ^ 



-3/4 



(7) 
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where uj is the Thompson cross-section and denotes N/10 x 
in cgs units. The large magnetic fields of the magnetar, decrease 
the effective cross-section (Herold 1979) but not sufficiently to 
make this radiation optically thin. With such a large optical 
depth the radiation forms a radiation-pairs plasma at a thermal 
equilibrium with an initial temperature of 
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300keV£i 6 %6 /2 Q 1/4 , (5) 
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where a is the Stephan-Boltzmann constant. This radiation- 
pairs plasma expands relativistically (Goodman 1986; Paczyn- 
ski 1986; Shemi & Piran 1990; Piran et al. 1993; Meszaros et al. 
1993; Katz 1996; Thompson & Duncan 1996) with T oc R and 
T oc R' 1 until at R± the pairs stop annihilating and their num- 
ber freezes. The number of remaining pairs, N±, is determined 



This kinetic energy is smaller by two orders of magnitude than 
the minimal energy required to produce the observed radio flux 
(Eq. 3). Therefore, the energy source of the radio afterglow 
cannot be the kinetic energy that remain in the pairs outflow. 

It is important to note that while strong magnetic fields, that 
might be dragged from the magnetar into the fireball, may influ- 
ence the interaction between the photons and the pairs they do 
not change the conclusion. Strong magnetic fields (B>10 13 G) 
would decrease the cross-section for photon-electron (positron) 
scattering (Herold 1979) resulting in a smaller T±. In addi- 
tion, a strong magnetic field would suppress the cross-section 
for pair annihilation into two photons, but it would also open a 
new channel of pair annihilation into a single photon (Wunner 
1979; Harding 1986). The latter becomes the dominant anni- 
hilation process and its cross-section is larger than the cross- 
section with no magnetic field that was used in Eq. 6. The 
overall result is a lower N±. Thus, strong magnetic field would 
only decrease the energy that remain in the pairs. The energy of 
the magnetic field itself may, however, contribute a significant 
component to the energy of the relativistic outflow. We address 
this contribution in section §4. 



3 This result corrects the one presented in Grimsrud & Wasserman (1998) that find no dependence of the final pairs energy on the initial radius. We show here that 
E± oc # 3 ^ 4 . The value of E± for Rq = 10 6 cm is similar in both works (see §5). 
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3.1. Interaction between the radiation and the circum burst 

medium 

The energy needed to power the afterglow is only a small 
fraction of the total prompt 7-ray energy. It is therefore worth- 
while exploring whether the interaction of the prompt radiation 
with the circum-flare medium can give rise to a relativistic out- 
flow with the required energy. The optical depth is given by 
GjnR, where n is the ambient density and R w 10 16 cm is the 
observed size of the radio emitting region. The energy acquired 
by the ambient electrons within R is therefore 

^151 = Gjn R = 5 x 1(TV? 16 (9) 
E 

Tapping 10" 3 of the energy requires an unreasonably large av- 
erage density of n > 10 5 cm" 3 . This is not expected around the 
magnetar given possible previous bursts. Moreover, giving this 
amount of energy to such a large mass would result in a sub- 
relativistic velocity (v -C 0.1c). 

An alternative mechanism can be pair enrichment by col- 
lisions of the outgoing radiation with photons that are back- 
scattered by the ambient medium (Thompson & Madau 2000; 
Beloborodov 2002; Meszaros et al. 2001). These pairs would 
in turn scatter more photons that would create more pairs and 
so forth. This process will take place up to a radius where each 
ambient electron scatters at least one photon: 



enrich : 



Ea T 
4ire~, 



1/2 



= 2.5 x 10 13 cm£4 6 /2 . 



(10) 



At smaller radii the number of pairs will grow exponentially 
until there will be about m p /m e pairs for each ambient medium 
electron. At this point there are enough pairs to accelerate the 
ambient medium to a relativistic velocity. Once this happen the 
temperature of the radiation in the rest frame of the accelerated 
medium drops significantly. Since the spectrum of the radiation 
is thermal there are no photons with an energy larger than m e c 2 
at this frame. The scattered photons stop creating additional 
pairs and the exponential process stops. Thus, up to R pa i rs the 
density of pairs is expected to be about m p /m e times larger than 
the external medium density. The fraction of energy acquired 
by this pair enriched region is 

£<■'<<■<■'■" ^ 7 rt„.,,,,,,- 10 '„„ (11) 



E m e 

Again, this enrichment is insufficient to tap 10~ 3 of the ini- 
tial energy unless the external medium average density is ~ 
10 5 cm" 3 . Even if one sets such large density around the neu- 
tron star, the high density region must be truncated shortly after 
3 x 10 13 cm, to allow for the observed mildly relativistic motion. 
Such a configuration seems to be too contrived. 

We conclude that the interaction between the prompt radia- 
tion and the external medium is unlikely to be the source of the 
afterglow energy. 

4. LOADED FIREBALLS 

4. 1 . Baryonic Load 

The processes that govern the evolution of a fireball loaded 
with protons are similar to those that govern a pure pairs- 
radiation fireball. However, the electrons that accompany the 
protons contribute to the opacity while the protons contribute 
to the inertia. Both effects can be taken into account by gener- 
alizing Eqs. 7-8. To do so we replace m e with the mean mass 
per particle: 

_ (m e /m p )N± + N p 



m = m„ 



N ± +N p 



(12) 



and we replace N± with the total electrons and positrons den- 
sity: 

N = N P +N±, (13) 

where N p = E/ (m p c 2 rj) is the number of protons and 77 = E /Mc 2 
characterize the baryonic load. The generalized equations are 
valid as long as the baryons load is small enough and the ra- 
diation escapes before most of its energy is converted to the 
baryonic kinetic energy. 

The behavior of m and Nb± as a function of the baryon load, 
77, depends on the ratio N p /N±. There are two critical values of 
this ratio: N p /N± = m e /m p and N p /N± = 1. The former marks 
equal mass for baryons and pairs, the latter marks equal contri- 
bution to the Thompson scattering. These values corresponds 
respectively to the following critical values of 77: 
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= 2.5xl0X 1 6 /4 <V 1 1/4 (15) 



An additional critical value of 77 is defined by the condition that 
the photons have transferred effectively all their energy to the 
baryons: 
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Eaj 



\4TTRoctm p c 2 
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= lOOElfR-ftl 4 . (16) 



For this critical 77 the photons decouple from the baryons just 
at the moment that the baryons stop being accelerated by the 
photons. 

Figure 1 illustrates the dependance of the energy that remain 
in the plasma, E Pi ± and its final Lorentz factor T ; , ± on 77. When 
rj\ <C 77 the baryons do not affect the evolution, m « m e and 
yV « N± and therefore the fireball evolves as a pure pairs-plasma 
fireball. The bulk Lorentz factor and energy are given by equa- 
tions 7 and 8. When 772 <C 77 <C 771 the baryons carry most of the 
inertia while the pairs are still responsible for most of the opac- 
ity. Therefore, m « m e rn/ri while N « N± and E Pt ± oc rj~ 3 ^ 4 
while its final Lorentz factor T Pt ± oc 77 1 / 4 . In this case the pairs 
become transparent for the radiation at R±. For 773 <C 77 <C 772 
the baryons provide both the inertia and the opacity. The final 
Lorentz factor is constant and E Pt ± oc 77"' . The radiation decou- 
ple from the electrons at radius larger than R±, but the radiation 
still carry most of the energy and its temperature is w 7b. For 
77 <C 773 the radiation decouples from the baryons only after it 
transferred most of its energy to the baryons. The baryons en- 
ergy is the total energy E and their Lorentz factor is 77. The 
radiation energy however is much smaller than E and its tem- 
perature is much smaller than 7b. We call this (77 < 773) a heavy 
load. 

The final Lorentz factor can be approximated as (see Fig 1): 



■p± 



680(7 ? /77 1 ) 1/4 £i 6 / X 1 6 

ioo4 /4 <v; /4 

77 



771 < 77 

V2<v<m 
V3<v<m 



(17) 



While the energy that remain in the ejecta is: 



Ep± 
E 



\.4 x 10-X 3 / 4 m<v 

4 x io-\r,/ V2 r 3/ %f V2<V<Vi 

mh Vi<v<V2 
1 77 < 773 



(18) 
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So far we have only considered protonic loading. The evolu- 
tion of a neutron rich fireball was explored in the past by several 
authors (Derishev et al. 1999; Beloborodov 2003; Rossi et al. 
2004; Vlahakis et al. 2003). The neutrons are coupled to the 
plasma only through collisions with protons. Initially the pro- 
tons drag the neutrons efficiently. A neutron that collides with 
a proton receives in such a collision at most lGeV. Therefore, 
the dragging takes place as long as each neutron collides with 
at least one proton during the time that the protons double their 
Lorentz factor. This happens when: 

Rn p a /T > 1 (19) 

where n p is the protons density and do = <J n pf3 r ei ~ 3 • 10~ 26 cm 2 
where a np is the neutron-proton cross section and j3 re i is the 
relative velocity between neutrons and protons. 

With (To about an order of magnitude lower than aj, the neu- 
trons always decouple from the protons before the radiation de- 
couples from the plasma. If the neutrons do not decouple be- 
fore the acceleration ends the load is essentially heavy and the 
photons transfer most of their energy to the baryons. The fi- 
nal Lorentz factor in this case is 77, corresponding to the total 
baryonic load. If the neutrons decouple before the acceleration 
ends then during the decoupling phase (when the neutrons be- 
gin to lag after the protons) the neutrons become relativistically 
hot and the inelastic n-n collisions convert neutrons to protons 
while producing pions. If initially N„/N p » 1, a significant 
fraction of the neutrons will be converted to protons and this 
ratio will become of order unity (Derishev et al. 1999) or some- 
what larger (Fuller et al. 2000) after decoupling. Thus the final 
protonic load, that carries most of the energy, is comparable to 
the initial total baryonic load. The final Lorentz factor and the 
energy of the plasma can be approximated by Eqs. 17 & 18 if 
77 is defined according to the total baryonic load. 

4.2. Heavy Loading 

For completeness we also consider here the energy and the 
temperature of the thermal radiation that escapes from a heavily 
loaded fireball, namely a fireball with 77 < 773. This situation is 
clearly inapplicable to this giant flare, but it may be relevant to 
GRBs, possibly as an explanation of precursors which are less 
energetic than the burst itself. In this case the Lorentz factor 
saturates at R n = Rot], and the fireball enters its coasting phase. 
The fireball is coasting without spreading (both in the local and 
the observer frame) until R spre ad = ctrj 1 . At larger radii the fire- 
ball spreads to a width of ~ Rjr\ 2 in the observer frame. The 
radiation in this case (77 < 773) decouples from the matter long 
after R v at: 



E(7j 
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(20) 



spread 



For Rot/ < R < mm{R p h, R spr ead), i-e during the coasting phase 
while the fireball is still opaque, V oc R 2 and therefore T oc R~ 2 I 3 
and E p h oc R~ 2 I 3 . If the fireball is still opaque at R sprea d then 
for Rspread < R < Rph the radiation evolves with T oc R~ l and 
E p h oc R~ x . Thus as far as the radiation is concerned there is 



another critical 77 for which R p h = R sp read- 



For 774 < 77 < 773: 
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For 1 < 77 < 774: 
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In GRBs the final Lorentz factor of the relativistic ejecta 
is >300, as indicated by opacity considerations (Lithwick & 
Sari 2001), and therefore ?7>300. On the other hand, the non- 
thermal spectrum of the prompt emission implies ?7<773 ~ 10 3 
for GRBs, making the range of allowed 77 very narrow. It 
implies also that the ratio between the energy in the ther- 
mal radiation that escapes the fireball and its kinetic energy 
is (77 3 / 77)- 8 / 3 >(1000/300)- 8 / 3 « 5%. Taking into account the 
efficiency in which this kinetic energy is converted into a 
non-thermal radiation, it indicates that there should be a non- 
negligible thermal component in almost any GRB. This result 
may be supported by the observations, as suggested by Ryde 
(2005). The idea that a thermal component in the prompt emis- 
sion is indeed the radiation that escapes from the fireball can be 
tested simply by comparing the energy and the temperature in 
this component: 



T ohs = lMeVE tht5l E^ /4 R Q6 

Where E t h is the energy in the thermal component and E is the 
total observed energy (note that we use here the typical values 
for GRBs). This relation should be tested within a given burst 
pulse by pulse. 

4.3. Electromagnetic load 

An alternative loading of the fireball is an electromagnetic 
load. Consider a magnetic field that is confined to the fireball 
and that is accelerated with it. The magnetic field carries an 
energy jE E m where E E m is the electromagnetic energy in the 
fireball's rest frame at the time that the pairs decouple from the 
radiation. The electromagnetic loading can be treated similarly 
to the baryonic loading where the magnetic field contribution 
to the inertial mass is: 

E 

m = m e + — 5, (25) 

VemN±c z 

and t]em =E/Eem is a measure of the electromagnetic loading. 
The total number of pairs, and therefore the opacity, is similar 
to those of a pure pairs-radiation fireball 4 : 

N = N±. (26) 

cross-section for Compton scattering is not affected. In our case where E » 10 46 



-1/2-1/4 
'-1 • 



(24) 



During the acceleration phase (R < R n ) the photons cool in 
their rest frame. The acceleration compensates for this cool- 
ing keeping the temperature and the overall radiation energy in 
the observer frame, E p h, constant. During the coasting phase, 
with no acceleration, the photons cool and lose energy in the 
observer frame as well. Since the ratio of the photon num- 
ber to the proton number, ~ m p c 2 rj/To 3> 1, is roughly con- 
stant at all times, the photons govern the cooling with an adi- 
abatic index of 4/3 (TV 1 ? 3 = const, where V is the volume). 

4 This is true as long as the magnetic field is much smaller than 10 13 G at R±, so the 
erg, R± fa 10 7 cm and if the relevant loading is t/em^IO this condition is satisfied. 
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FIG. 1 . — The baryons' energy (thick line, left y-axis) and the baryonsS final Lorentz factor (thin line, right y-axis) as functions of r\ (lower x-axis). For clarity the 
total baryons mass is depicted on the upper x-axis. Also marked (short dashed horizontal lines) are the lower limit on the energy inferred from the radio afterglow 
and the energy that remains in the pairs in a pure radiation-pairs fireball. The vertical dashed lines mark the allowed loading in order to produce a thermal flare at 
the observed temperature. The parameters in this plot are E = 3.5 X 10 46 erg, t = O.lsec and Rq = 10 6 cm. 



Therefore, there are only two critical values of tjem- Vema = f]\ 
and rjEM.i = Tq/T±. When 771 <C t\em the fireball evolves as a 
pure pairs-plasma fireball. When t)em,2 -C t]em <C t]em,\ the ra- 
diation still carries most of the energy and its observed temper- 
ature is To. The total energy that is left in the magnetized pairs 
after they decouple from the radiation is Eem.± <x rf^l A . Finally, 
when rjEM <SC rjEM.2 the radiation decouples from the magne- 
tized pairs only after it has transferred to them most of the en- 
ergy. In this case Eem.± = E while The radiation energy is E p h = 

E(7]EM,2/riT &/3 and its temperature T obs = T (r]EM,2/rjT S/3 - Not e 
that when the electromagnetic loading is significant (rjEM <Vi) 
the Lorentz factor of the pairs at the time that they decouple 
from the radiation is only a lower bound on their final Lorentz 
factor. The reason is that energy transfer between the electro- 
magnetic field and the pairs can still take place at larger radii. 

5. COMPARISON WITH PREVIOUS WORKS 

The evolution of pure and loaded fireballs was explored 
by several authors in the past. Paczynski (1986) and Good- 
man (1986) considered only the properties of the radiation that 
emerges from a pure fireball. Shemi & Piran (1990), Piran et al. 
(1993) and Meszaros et al. (1993) considered the evolution of 
a loaded fireball. These papers miscalculated the radius that 
the baryons decouple from the radiation, resulting in an over- 
estimate of 773 and of the final Lorentz factor and the ejecta en- 
ergy, for 77 > 773. In the context of GRBs our result (773 ~ 10 3 ) 
is significantly lower than in those papers, and combined with 
the lower limits on the Lorentz factor, 77 500 is narrowly con- 
strained. More recently, Grimsrud & Wasserman (1998) carried 
out a detailed numerical and analytical calculation of the final 
Lorentz factor and energy of pure and loaded fireballs. Our 
much simpler estimates generally agree with their results with 
the exception of the final energy that remain in the pairs, E± and 



their opacity for the radiation once the Lorentz factor saturates 
T±. While Grimsrud & Wasserman (1998) find that both quanti- 

-3/4 

ties do not depend on the initial radius, we find that E± oc R 

-3/4 

and that t± oc R as well. However, the value that these au- 
thors obtain for Rq = 10 6 cm are similar to the values that we 
present here. 

6. DISCUSSION 

We presented here a simple, yet comprehensive, derivation 
of the evolution of pure and loaded fireballs. Our results cor- 
rect previous works (see §5), enabling us to compare the energy 
that remains in the plasma once the radiation escapes to the ob- 
servations of the radio afterglow of the giant flare from SGR 
1806-20. 

The relativistic ejecta that produced this radio afterglow con- 
tained at least 0.1% of the energy emitted in 7-rays. The ther- 
mal spectrum of the flare and its temperature indicate that the 
flare resulted from an energy deposition near the neutron star 
surface. Such energy deposition must create an opaque acceler- 
ating fireball. An elegant explanation for the afterglow energy 
source could have been the inevitable energy of the remaining 
pairs after they decouple from the radiation. However, we find 
that this energy is short by at least two orders of magnitude, 
excluding this possibility. We considered also the possibility 
that the energy of the afterglow is obtained by the interaction of 
the outgoing radiation with the external medium. We find this 
scenario to be unlikely since it requires a rather high external 
density with a contrived profile. 

We conclude, therefore, that if the relativistic ejecta that pro- 
duces the afterglow is directly related to the prompt flare (rather 
than e.g. to the following confined fireball), then the fireball 
must have been loaded with either baryons or Magnetic field in 
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order to enhance the energy that remains in this ejecta. This 
loading should be however fine tuned in order to obtain the 
right amount of energy. Even if the afterglow energy is com- 
parable to that of the ejecta (we have only lower limit on this 
energy) the loading cannot be too high. If r\ < 773 the energy in 
the flare drops significantly and so does the temperature. The 
agreement of the observed temperature and the black body tem- 
perature estimated in Eq. 5 suggest that this is not the case. For 
the December 21 th giant flare the range of the allowed baryonic 
loading is 100<?7< 10 5 , which corresponds to a baryons mass 
of 10 20 gr^^l0 23 gr(seefig. 1). 

Interestingly similar lower limit on the ratio of the afterglow 
energy to the flare energy was observed in the August giant 
flare from SGR 1900+14 (Frail et al. 1999), although this flare 
was hundred times dimmer than the December 21 th giant flare. 
With the caveats of small number statistics, and both afterglow 
estimates being lower limit, this similarity suggests a common 
origin for the two afterglows and a linear relation between the 
observed afterglow energy and the flare energy. The similar- 
ity between the fractional energy in the afterglows of the two 
events, is not accounted for in our model. The most natural 
energy source, the pairs energy in a baryon free fireball, con- 
tain a constant fraction of the fireball energy (see Eq. 8) but 
this fraction is far too short. If this similarity is confirmed by 



future events it may indicate that some process regulates the 
baryonic load. For example it is possible that a constant num- 
ber of baryons per unit energy are torn apart from the surface of 
the neutron star and are mixed with the fireball. Alternatively 
it is possible that a fixed fraction of the energy is ejected as a 
Poynting flux. A different solution might be that the fireball 
is heavily loaded aspherically, such that only a small portion 
of the 4-n solid angle is loaded, while the rest of the fireball is 
pure. In this case E P: ± /E is simply the ratio of loaded portion 
solid angle to the pure portion solid angle, and this ratio might 
be similar between different flares. It is, of course, also possible 
that the afterglow is not related directly to the prompt flare and 
that it is created by an independent mechanism and its energy 
is not extracted of the initial fireball. Once more the similarity 
between the ratios of afterglow to flare energies in both cases 
poses a puzzle for this last explanation as well. 
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